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ABSTRACT 

We describe deep, new, wide-field radio continuum observations of the Great Observatories 
Origins Deep Survey - Nortlr (GOODS-N) field. The resulting map has a synthesized beamsize 
of ^1.7" and an r.m.s. noise level of ■~3.9/^Jy beam~^ near its center and ^8/iJybeam^^ at 15' 
from phase center. We have cataloged 1,230 discrete radio emitters, within a 40' x 40'region, 
ab ove a 5-a detectio n thre shold of '^20/xJy at the field center. New techniques, pioneered 
bv lOwen fc MorrisonI ( 20081) . have enabled us to achieve a dynamic range of 6800:1 in a field 
that has significantly strong confusing sources. We compare the 1.4-GHz (20-cm) source counts 
with those from other published radio surveys. Our differential counts are nearly Euclidean be- 
low 100 /xJy with a medi an source diameter of ^1.2". This adds to the evidence presented by 
Owen &: MorrisonI (j2008() that the natural confusion limit may lie near 1 /iJy If the Euclidean 
slope of the counts continues down to the natural confusion limit as an extrapolation of our log N 
- log S, this indicates that the cutoff must be fairly sharp below 1 /iJy else the cosmic microwave 
background temperature would increase above 2.7K at 1.4 GHz. 



Subject headings: galaxies: evolution - 
galaxies 

Introduction 



- galaxies: galaxies: starburst: galaxies: AGN — radio continuum: 



The GOO DS-N field (|Dickinson et al.l 120031: 



Giavalisco et a l. 2004i) covers ss 160a rcmin^ cen 



tered on the Hubble Deep Field North (jWilliams et al 
19961 ) and is unrivaled in terms of its ancillary 
data. These include extremely deep Chandra, 
Hubble Space Telescope and Spitzer observations, 
deep UBVRIJHK ground-based imaging and 
^3,500 spectroscopic redshifts from 8-10-m tele- 
scopes. Previous radio observations of this region, 
however, fell short of complementing this unique 
dataset. 

Radio emission is a relatively unbiased tracer 
of star formation and can probe heavily obscured 
active galactic nuclei (AGN) - objects that are 
missed by even the deepest X-ray surveys. Ra- 
dio observations thus allow us to fully exploit the 
wealth of data taken at X-ray-through-millimeter 
wavelengths, providing a unique extinction- free 
probe of galaxy growth and evolution through the 
detection of starbursts and AGN. 



The recent imaging of lOwen fc MorrisonI (|2008l ) 
(OM08) (cr = 2.7AiJybeam-i at 1.4GHz) have 
shown that the techniques exist to make radio im- 
ages that approach the theoretical noise limit. To 
this end, we have obtained new, deep radio imag- 
ing of the GOODS-N field. 

While GOODS-N was selected to be free from 
bright sources at optical wavelengths, the field 
contains several very bright radio sources which 
place severe limitations on the dynamic range that 
can be obtained and hence the ultimate sensitiv- 
ity of the radio map. Before new techniques were 
developed to deal with this issue, only mo derately- 
deep radio imaging was possible ([Richa rds 2000"). 

The earliest VLA data in the GOODS-N 
fie ld were re processed using ne w techniques b y 
.Biggs fc Ivison (.2006 ') and , Morrison et al.l (|2008l ). 



achieving a noise level of 5.7-5.8 /it Jy beam ^ - 
a 23-25% impr ovement on the original map by 
Richards! (l2000l) . and close to the theoretical noise 



limit when one considers the increase in system 
temperature at the low elevations permitted dur- 
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Table 1 

GOODS-NoRTH VLA observing log 



Observing dates 


Integration time (h) 


NRAO project code 


Configuration 


Nov-Dec 1996 


42 


AR368 


A 


Feb 2005 


28 


AM825 


B 


Aug 2005 


7 


AM825 


C 


Dec 2005 


2 


AM825 


D 


Feb-Apr 2006 


86 


AM857 


A 



ing the original observations. 

While the reduction of Biggs fc Ivison ( 20061 ) 
provided improved access to the ^Jy radio pop- 
ulation, even deeper radio imaging is required to 
properly complement the extremely deep GOODS 
Spitzer mid-infrared data, as well as forthcom- 
ing deep observations at far-infrared and submil- 
limeter wavelengths from Herschel, SCUBA-2, the 
Large Millimeter Telescope, and other facilities. 
To this end, we have added 123 hr to the ex- 
isting data. The reduced full resolution image 
(beam=1.7") and its rms map are available on- 
lin4!|. 

The paper is laid out as follows: in §2 we de- 
scribe the observations and the reduction of the 
data. In §3 we discuss the cataloging of radio emit- 
ters. §4 contains the results and some discussion 
of the catalog. We present our conclusions in §5. 

2. Observations, Reductions, and Cata- 
loging 



In 1996 November, Richardj (|200 d) observed a 
region centered at 12:36:49.4, -^62:12:58 (J2000) 
for a total of 50 hr at 1.4 GHz using the National 
Radio Astronomy Observatory's (NRAO's) VL^Il 
in its A configurat ion. Of this, onl y 42 hr was con- 
sidered usable by iRichardd ( 2000f ). Adopting the 
same position and frequency, we obtained 28 hr of 
data in the VLA's B configuration in February- 



^http:/www. ifa.hawaii.edu/~morrison/GOODSN and at 
the NASA/IPAC Infrared Science Archive (IRSA) 
^ttp://irsa. ipac.caltech.edu/ as an ancillary data product 
associated with the GOODS Spitzer Legacy survey 

*The NRAO is operated by Associated Universities, Inc., 
under a cooperative agreement with the National Science 
Foundation. 



April 2005, 7hr in C configuration in August 2005, 
and 2hr in D configuration in December 2005, 
and 86 hr in A configuration in 2006 February- 
April (see Table [ij - for a useful combined to- 
tal of 165 hr. Observations were done at night to 
avoid solar interference. We followed the 1:4 scal- 
ing of visibility data between the arrays described 
by OM08. This empirically derived scaling rela- 
tion provides for more uniform weighting of uv 
data. 

In most regards the new observations were 
taken using the same parameters as those used in 
1996. However, the integration time was changed 
from 3.33 s to 5 s because of difficulties experienced 
by the correlato r with the shorter integration time 
(lRichardsll207)0h . The data were all obtained using 
spectral-line mode 4, which yields 7 x 3.125-MHz 
channels in each of two intermediate frequencies 
(IFs), centered at 1,365 and 1,435 MHz, in each of 
two circular polarizations. The channel width and 
integration time were compromises chosen to max- 
imize reliability and sensitivity while minimizing 
bandwidth (radial) and time (tangential) smearing 
effects, respectively. The upcoming EVLA cor- 
relator, WIDAR, will offer much shorter integra- 
tion times, narrower channels, and greater overall 
bandwidth. 

2.1. Calibration and editing 

AXVS was used to reduce and analyze all the 
radio data. The first step was the calculation and 
corrections of the spectral bandpass shape. This 
was done in the following manner using the task 
BPASS. The bright point-source phase calibrator 
was SPLIT from the raw database and phase self- 
calibration was applied. The self-calibrated data 
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Cont peak flux = 6.7631 E-05 JY/BEAM 
Levs = 1.000E-06 * (4, 4.500, 5, 7, 10, 20) 



Fig. 1. — Contours of constant r.m.s. noise after correction for the shape of the primary beam. The GOODS-N 
area covered by HST/ ACS is outhned near the field center, with an average sensitivity of '^3.9^Jybeam~^. 
At a radius of ~15' (roughly the half-power point of the primary beam) the r.m.s. sensitivity drops to 
~7 /iJybeam"^. The strong inward defiection of the 7-/iJy and 10-/iJy contours is due to the 265-mJy radio 
source discussed in ij4.4.1l 



were then used to calculate a bandpass correction calibrated. The uv data for the target field were 

and flatten the spectral response across the band split from the database and clipped using the 

for the uncalibrated multi-source database. Stan- AXVS task clip at a level well above the total 

dard flux density calibration was applied next, flux density found in the fleld to remove any inter- 

using the Baars flux-density scale with 3C286 ference prior to the self-calibration process. Only 




minor interference was encountered during the ob- 
servations. 



based weights for each 5-s integration were also 
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2.1.1. UVFIX - Calculating w, u, and w 

In verifying the astrometry of the A-array data 
from 1996 and then again in 2006, we found a ro- 
tation between the two astrometric frames. The 
rotational offset was about 1" at a radial distance 
of 20' from the phase center, and is hkely the cause 
of the optical-radio offsets reported reg ularly over 



the last decade (e.s.. Ilvison et al.ll2002[ ). The off- 



set problem was traced back to the VLA on-line 
system (the so-called MODCOMPS) and its cal- 
culated u,w, and w terms. This problem was cor- 
rected by using the AIPS task UVFIX. This task 
computes values of u,v, and w using the time in- 
dex and baseline information recorded during the 
observations. The AIVS task UVFIX was used 
on all the uv datasets after the initial bandpass 
and flux calibration. 

2.2. Imaging and self-calibration 

Bright sources and their sidelobe patterns affect 
the noise structure and the dynamic range of the 
final reduced maps if they are not cleaned. Such 
sources were located using D-configuration data. 
We imaged to a radius of 8 degrees using a pixel 
size of 30". The GOODS-N field was chosen to 
be free from bright objects at most wavelengths, 
but not the radio waveband. As a result, bright 
radio sources limit the dynamic range. We imaged 
sources as far as 1.5 degrees from the phase center 
- a total of 56 facets, selected to cover distant, 
bright sources, plus another five facets within the 
primary beam to help with the PEELR process 
which is fully discussed in OM08. 

Each facet is constructed from the Fourier 
transform of the data, phase shifted to its tan- 
gent point on the celestial sphere. In all, 98 facets 
were used. This technique is known as polyhedral 
imaging (als o known as faceting) and is discussed 
at length bv lCornwell fc PerlevI (|l992l ). 

One night of A-configuration data was used to 
lay out 37 facets to cover the primary beam at 
1.4 GHz. Each of the 37 facets comprised 1024^ 
pixels, each 0.5" pixels. The other 56 facets were 
512^-pixels in size, again with 0.5^"pixels, cen- 
tered on the bright, outer sources. 

Initial cleaning of the fields with IMAGR re- 
vealed faint sources within the 98 facets and tight 
clean boxes were placed around each one, thus lim- 
iting cleaning to real features. Maps were cleaned 



down to the l-a noise level to create clean compo- 
nents for use in the self-calibration process. Clean 
boxes are a vital aspect of this process. They avoid 
the possibility of cleaning into the noise, which can 
change the noise structure of a map and result in 
so-called 'clean bias'. Failure to use clean boxes 
can artificially reduce the noise. 

This one night of uv data were first self- 
calibrated in phase only, via the AXVS task 
CALIB, and a new map was made. Self- 
calibration in amplitude and phase followed. The 
clean components from this one night were used 
to construct the fiducial model for the remaining 
13 nights of new A-configuration data (AM857 - 
NRAO project number) as well as the seven nights 
of 1996 (AR368) A-configuration data. For the 
four nights of new B-configuration data (AM825), 
one of the B-configuration datasets were used to 
provide the clean model for self-calibration. The 
C- and D-configuration data (AM825) were boot- 
strapped to the clean model in a similar fashion. 

Editing of the uv data was done before and after 
the self-calibration process using the ALPS task 
TVFLG. After each episode of self-calibration, 
the clean component model was removed from the 
uv database using UVSUB. UVPLT was then 
used to determine the residual amplitude distribu- 
tion. Any high residual points were deleted with 
CLIP, except those on the shortest baselines. 

The 165 hr of data resulted in a very large 
database so we reduced the data volume by ex- 
ploiting the fact that the A- and B-configuration 
(programs AM825 and AM857) observations were 
taken at the same hour angles (±3.5 hr of transit) 
over many nights. The smaller hour angle inter- 
val helped mitigate the higher system temperature 
that occurred at lower dish elevations. This is due 
to the fact that the antenna feeds see the ground at 
low elevations. The AR368 observing tracks cov- 
ered ±7.5 hr of transit yielding extended periods of 
low-elevation observations. The antenna weights 
were calibrated as part of the standard calibration 
process minimizing the impact of the more noisy 
low-elevation data. Finally, STUFFR is used to 
average data within a small range of u,v (so there 
is no effect on the image), which was weighted 
properly. In this way the number of data points 
used for the imaging were reduced by an order-of- 
magnitude. This is fully described in OM08. 

The five separate datasets were combined into 
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one large uv dataset using DBCON. This was 
then divided up by hour angle, by Stokes parame- 
ter (LL and RR, separately) and by IF, thus yield- 
ing eight separate datasets. This last step mini- 
mizes the problem of 'beam squint' wherein the 
two circular polarizations have different pointing 
centers due to the slight off-axis location of the 
feeds. This, combined with the alt-az telescope 
mount, causes the effective gain for an off-axis 
source to vary as a function of time. Moreover, the 
two IFs are at different frequencies and so yield a 
different primary beam shapes hence correspond- 
ingly different gains for off-axis sources. 

After self-calibration, some artifacts associated 
with two bright sources (within the 37 facets 
that map out the primary beam) remained. A 
local self-calibration technique in AXPS called, 

PEELR WAS USED TO REMOVE THESE SIDE- 
LOBES. The method is described as fol- 
lows. 

In brief, PEELR subtracts the clean 

COMPONENT (CC) model FOR ALL FACETS 
FROM THE UV DATA, EXCEPT FOR THE FACET 
CONTAINING THE SOURCE TO BE ' PEELED'. THEN 
PEELR SELF-CALIBRATES USING THIS INFOR- 
MATION (flux) in this facet, writes a new 

CALIBRATION TABLE. NEXT, PEELR REMOVES 
THE CCS FOR THE CHOSEN FIELD, THEREBY 
REMOVING THE BRIGHT SOURCE. THE SPECIAL 
CALIBRATION FOR THIS ONE FIELD AIDS IN THE 
SIDELOBE CORRECTION PROCESS. PEELR THEN 
GOES BACK TO THE ORIGINAL CALIBRATION TA- 
BLES AND ADDS THE CCS BACK TO THE UV 

DATASET. Having the uv data in eight sepa- 
rate FILES ALLOWS FOR A MORE ACCURATE CC 
MODEL, ENABLING MORE ACCURATE REMOVAL 
OF THESE BRIGHT SOURCES. 

Full details on this important process can be 
found in OM08. 

2.2.1. Final signal and noise images 

The final maps were constructed by combin- 
ing the 37 central facets made with each of the 
eight UV datasets, weighting them by 1/cr^, us- 
ing FLATN. The r.m.s., before correcting for the 
shape of the primary beam, is ~3.9 /iJy over a re- 
gion of 100^ pixels. The synthesized beam size 
is 1.7" X 1.6" with a position angle of —5 degrees. 
The final step was to run MVS task MWFLT 



which low-pass filters the image with a 101^-pixel 
kernel to yield a more uniform background. 

The AXVS task, RMSD, was used to construct 
a noise image, calculating a histogram based r.m.s. 
for each pixel using the surrounding pixels within 
a radius of 100 pixels. A multi-iteration process 
rejects pixels outside of the ±3(T range, thereby of- 
fering a more robust estimate of the noise. Fig. [1] 
shows contours of constant noise over the central 
region, after correcting for the primary beam re- 
sponse. 

3. Cataloging 

3.1. Angular size effects 

The angular size of discrete sources in the im- 
age are broadened by three effects: (1) the finite 
bandwidth of each channel where 



BWS '-^ — X (?•), (1) 

where /S.v is the fractional bandwidth (i.e. 3.125 MHz) 
and r is the distance from the phase center; (2) 
is the finite data sampling rate (estimated at a 
few percent); (3) and the true angular size of the 
source. Thus, in order to detect sources above 
4.5-sigma, images at 3"and 6"are also needed to 
detect extended sources-whether truly extended 
or experimentally extended. OM08 showed that 
3" and 6" maps are a useful complement to the 
full-resolution images to recover the full flux of ex- 
tended sources. The 3" and 6" images were made 
using the MVS task, CONVL. 

3.2. Source extraction and cataloging 

The AZVS task, SAD, ('search and destroy') 
was used to generate an initial source catalog. 
For each resolution, SAD was used in 'signal-to- 
noise ratio (SNR) mode' to search for peaks more 
than 4.5 X the local noise and to correct for radial 
smearing and primary beam attenuation. This 
mode uses the noise maps created by RMSD. The 
final catalog was clipped at >5f7. 

Residual maps created by SAD were searched 
by eye to find sources missed by the automatic 
procedure. We made SNR maps from the residual 
maps, searching for peaks above four. This lower 
SNR limit was picked because BWS reduces the 
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Table 2 

Sample table of radio sources 



Number 


R.A. (.12000) 


Dec. 


SNR 


Sp 


s. 




Size 




Upper 


Beam 




hms it e(s) 


d ' "± e(") 






(A'Jy) 


Maj 


Mill 


P.A. 


(") 


(") 


5 


12 34 3.51 0.061 


62 14 20.7 0.03 


27.6 


191. 6± 6.9 


795. 9± 44.9 


2.6 


0.7 


95 


0.0 


1.7 


9 


12 34 8.79 0.164 


62 9 20.5 0.10 


8.0 


51. 6± 6.5 


142. 4± 19.5 


0.0 


0.0 





3.8 


1.7 


12 


12 34 9.23 0.098 


61 56 45.5 0.09 


44.3 


397. 0± 8.9 


2767. 1± 171.5 


7.2 


0.0 


129 


0.0 


6.0 


13 


12 34 9.87 0.238 


62 3 57.6 0.12 


7.3 


57. 3± 7.9 


208. 6± 30.4 


2.8 


0.0 


89 


0.0 


1.7 


14 


12 34 10.61 0.085 


62 6 15.8 0.05 


18.6 


132. 1± 7.1 


564. 4± 43.5 


2.4 


0.8 


67 


0.0 


1.7 


18 


12 34 10.78 0.083 


62 2 54.8 0.05 


19.2 


173. 4± 9.0 


848. 0± 63.7 


2.6 


0.5 


63 


0.0 


1.7 


22 


12 34 11.16 0.029 


62 16 17.3 0.02 


239.1 


1674. 0± 7.0 


7741. 0± 236.7 


8.5 


1.6 


126 


0.0 


6.0 


23 


12 34 11.74 0.003 


61 58 32.5 0.00 


438.5 


4385. 0± 10.0 


27684. 0± 834.7 


1.7 


0.6 


52 


0.0 


1.7 


35 


12 34 16.09 0.071 


62 7 22.4 0.05 


22.8 


127. 5± 5.6 


346. 2± 18.2 


0.0 


0.0 





3.7 


3.0 


36 


12 34 16.17 0.157 


62 25 58.8 0.13 


8.6 


66. 1± 7.7 


284. 2± 34.6 


0.0 


0.0 





3.6 


1.7 



peak flux of a source. Next, the properties of these 
sources were measured using AIVS task, JMFIT. 
Sources with SNR greater than 4.5 were retained. 
Any such sources were added to the catalog for the 
appropriate resolution. 

SAD sources with SNR < 5.5 were re-measured 
by hand using JMFIT - the same task used by 
SAD to determine the source properties. At this 
low flux limit SAd's automatic detection routine 
does not always choose the optimal area for anal- 
ysis. Those pixels within a defined radius (100 
pixels) were used to calculate the local noise for 
each source. 

Following OM08, resolved sources were classi- 
fied on the basis of the best-Gaussian-fit major 
axis. If the lower limit for this parameter was 
greater than zero, the source is classed as resolved 
and the integrated flux was used for the total flux. 
If lower limit for the major axis was equal to zero 
(i.e., unresolved) then the peak flux was adopted 
as the total flux for the source. 

There are several extended sources within the 
central field, as shown in Fig. H TVSTAT was 
used to measure the total flux density of extended 
sources. The SNRs of extended sources were de- 
termined as follows: IMEAN was used to deter- 
mine the brightest peak and its position was then 
examined in the noise map; the ratio of these val- 
ues was adopted as SNR. 

Results were then collated following the pre- 
scription given in OM08. An additional noise term 
of 3 % was included for each source to account for 
calibration errors and uncertainties in the primary 
beam correction. 



4. Results and discussion 

4.1. Radio catalog 

A sample of the radio catalog is given in Ta- 
ble [2] The complete catalog is available electroni- 
cally. Column 1 is the source number: source num- 
bers below 10,000 relate to those sources found by 
SAD while those above 10,000 were those found 
and investigated by hand. Columns 2 and 3 con- 
tain the right ascension (R.A.) and the declination 
(Dec.) in J2000, with the positional error. Col- 
umn 4 represents the SNR using the ratio of the 
observed peak flux density and the local noise. In 
columns 5 we list the peak flux density and its 
error, in /ijybeani^^. Columns 6 contain the to- 
tal, primary-beam-corrected flux density and un- 
certainty, in ^Jy. The best-flt deconvolved size (in 
arcsec) is given in column 7-9. Using our criteria 
for resolved sources, if a resolved 2-dimensional 
Gaussian size was the best flt then we give the 
resulting major and minor FWHM size (in arcsec) 
along with the position angle (P.A.). If this was 
not the case then in column 10 we report the major 
axis upper limit, as estimated by JMFIT or SAD. 
If the source is resolved then we report 0.0" for the 
upper limit size. As noted above, very extended 
sources, as shown in Fig.lSllTTl had their positions, 
sizes, fluxes and errors measured interactively us- 
ing IMVAL and TVSTAT. The final column re- 
ports the resolution of the map that was used for 
the measurements for that particular source. 

4.2. Angular size distribution 

As a check on the radio angular size distribu- 
tion in the GOODS-N field, we have repeated the 
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Table 3 
Source size summary 



Radius 


Minimum 


Maximum 


Sources 


Upper 


Mean 


Error 


Median 


Resolved 


(arcmin) 


(^Jy) 






limit 


(arcsec) 


(arcsec) 


(arcsec) 


(%) 


5 


20 


30 


32 


37 


1.2 


0.1 


1.16 


46 


5 


30 


100 


49 


24 


1.8 


0.2 


1.15 


67 


10 


100 


300 


34 


26 


2.3 


0.4 


1.09 


57 


20 


300 


1000 


48 


11 


3.4 


0.6 


1.57 


81 



source size analysis of OM08. The focus of this 
analysis is to estimate the angular size of the faint 
radio population. As noted in column 1 of Table[3l 
we have analyzed the more sensitive inner region 
(<5') of the radio map for fainter sources while for 
the brighter sources we expanded our search out- 
ward (to 20') to improve the statistics. Columns 
2 and 3 give the flux range and column 4 lists the 
sources within that size range. Column 5 shows 
the number of upper limits used with the mea- 
sured source size distribution in the Kaplan-Meier 
analysif]^ to derive a robust mean size (column 6) 
and its associated error (column 7). The median 
size is reported in column 8 and the percentage of 
resolved sources used in this analysis is listed in 
column 9. 

The res ults in Table [51 are consistent with both 
OM08 and lMuxlow et all (2005 ). adding to the ev- 
idence that the size distribution of the radio popu- 
lation does not continue to decrease below 100 /xJy 
and that the median size of such faint emitters is 
around 1". 

4.3. Log iV Log S 

We have determined the differential normalized 
sources counts from our radio catalog using the 
code and analysis techniques described in detail 
by OM08. The resuhs are tabulated in Table E] 
The relative point-source sensitivity as a function 
of the distance from the field center for the three 
resolutions can be seen in figure 4 of OM08. Their 
figure 5 also shows the relative sensitivity as a 
function of source size. To summarize, multiple 



effects can make the source list incomplete, even 
at rather large SNR. The sensitivity for resolved 
sources is thus a function of both beam size and 
distance from the phase center. These effects in- 
clude (i) primary beam attenuation, (ii) reduction 
of peak flux due to BWS, (iii) finite time-average 
smearing and (iv) resolution bias. Sources with in- 
tegrated fluxes near the survey limit, if extended, 
will have peak sources below the limit. 

As seen in Table |3l the median source size is 
similar to that of the synthesized beam (~1.7"). 
To correct for the effects of this, we adopt the 
OM08 size distribution, which is based on the 
properties of sources at higher SNR. 

4.4. Comparison to previous surveys 

4.4-. 1. Dynamic range 



^IRAF task kmestim ate and algorithm by 
iFeigelson fc Nelsod lligSSl ) 



As described by Biggs fc Ivison ( 2006[ ). and 
seen in their flgure 2, the brightest GOODS-N ra- 
dio source within the primary beam of the VLA is 
near 12:34:52, -|-62:02:35. With a primary-beam- 
corrected total fiux of 263 mJy and a peak flux of 
~86mJy beam"^, it is bright sources like this one 
(see Fig. [10] upper right) that can limit the dy- 
namic range in an image. This is because there is a 
limit to how well we can clean these bright sources 
- a limit determined by how well we can character- 
ize the bandpass, the pointing, the antenna-based 
gains, and many other parameters. As a result of 
the imperfect cleaning, residual fiux spills into the 
image from the sidelobes of bright sources. Peak 
flux density, in Jybeam"^, is therefore a key quan- 
tity when determining the suitability of a held for 
deep radio imaging, where the 'beam~^' indicates 
a dependence on resolution. 
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Table 4 

Differential normalized source counts for GOODS-N 



5i 


5h 


(A^Jy) 


No. 

# 


(Jyi-^sr-i) 


28.0 


40.0 


34.0 


185 


7.37±0.70 


40.0 


55.0 


46.7 


187 


7.45±0.89 


55.0 


90.0 


71.1 


274 


6.67±0.70 


90.0 


150.0 


113.1 


150 


7.39±1.31 


150.0 


300.0 


199.8 


107 


5.29±1.11 


300.0 


1,500.0 


557.3 


65 


6.02±1.79 


1,500.0 


60,000.0 


5,260.5 


17 


5.56±1.68 



The reduction by iRichardsl (l2000l) yielded an 
image with a dynamic range of ~ 2, 800 : 1, with 
a = 7.5 fiJy. This was due in part to the fact that 
self-calibration was not apphed to the uv data. 
Self-cah bration and newly de veloped techniques 
allowed Biggs fc IvisonI ( 2006h to create a more 
sensitive map, with cr = 5.8 fxJy and a dynamic 
range of ~ 4,400 : 1. Our reduction, which in- 
cludes a large quantity of new data (in all VLA 
configurations) , reaches a = 3.9/iJy beam~^ and 
a dynamic range of ~ 6, 800 : 1. 

Th e ~140-hr VLA observation bv lOwen fc MorrisonI 
(|2008l ) - in 1046+59, a field deliberately chosen to 
be free of strong radio sources - has a = 2.7 fi3y, 
which is close to the theoretical noise limit, yet 
the dynamic range is only ^ 2, 500 : 1 because of 
the paucity of bright radio emitters. 

44.2. Log N-log S 

Fig.[2]shows our differential source counts below 
1 mjy, compared to other deep, 1.4- GHz contin- 
uum surveys. Our source counts agree with those 
for the 1046-1-59 field where the same methodology 
was employed for cataloging and for determining 
log iV-log S: use of multi-resolution maps to build 
a master source catalog, and the adoption of the 
same source size distribution. 

Above 100 /iJy o ur new survey agrees w ith the 
counts presented bv lBiggs fc Ivison (l2006l) . Below 
100 /iJy, the counts in the new surv ey flatten out 
while those of Biggs fc IvisonI ( 2006 ) began to turn 
down. This is a well-known issue in deep radio 
surveys. The difference can be related to a number 



of complications: 

1. We believe that one of the main problems 
arises from the Gaussian fit not always being 
the best approach to estimate source prop- 
erties. Indeed, the intrinsic shape of the 
sources gets badly affected due to bandwidth 
and time delay smearing. We have found ev- 
idence that flux densities are underestimated 
systematically after sources have been con- 
volved to lower resolutions - flux densities 
are larger in the convolved maps. Also, a 
particular uncertainty comes from the as- 
sumption that in unresolved sources we con- 
sider the peak flux density as an estimate 
of the total flux density. These problems 
clearly suggest the necessity for a more ro- 
bust source extraction method in radio ob- 
servations. Another, smaller uncertainties 
coming from the 'clean bias' and 'flux boost- 
ing' are not expected to be important due to 
the tight clean boxes we used, and the high 
5-CT threshold in the source extraction, re- 
spectively. 

2. A knowledge for the intrinsic radio source 
size distribution is essential when estimating 
the number of sources miss e d by o ur imaging 
approach. IMuxIow et al.l (|2005[ ) use high- 
resolution MERLIN and A-conflguration 
VLA data to tackle this problem, but they 
miss a non negligible number of low-surface- 
brightness sources, as is indicated clearly by 
cross-matching the catalog with lower res- 
olution observations using the WSRT. We 
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believe the combination of VLA data in its 
A, B, C, and D configurations followed by 
a source extraction performed in convolved 
images is essential for detecting these low 
surface brightness sources (>4"). A simple 
source extraction using SAD in the original 
non-convolved image misses a considerable 
number of sources as we show in this work. 
We find that our survey is in good agree- 
ment with OM08, so given by the similar 
treatment, we have adopted their source size 
distribution for number count estimations. 

3. Finally, the source extraction based on a 
SNR criterion highly depends on reliable es- 
timates of noise as a function of position 
in the map and the efficiency for detection. 
Monte-Carlo simulations for the source ex- 
traction method (e.g. inserting randomly 
placed sources in the map and then doing 
a random search to extract them) are typi- 
cally used for these purposes. Nevertheless, 
given our detailed source extraction and our 
constant examination of residual maps, we 
believe that our number counts do not suffer 
significantly from incompleteness down to a 
threshold of 5a. This is why we prefer to 
perform theoretical corrections, as explained 
in OM08, rather th an a simple Mont e-Carlo 
approach as was in llbar et al. I (l2009l) . 



4- 4- 3- Log N -log S Euclidean nature below 
lOOfiJy 

It is expected that the Cosmic microwave back- 
ground (CMB) temperature would increase above 
2.7K at 1.4 GHz if the counts do not turn over at 
some flux density below where we have reached. 



(e.g. iGervasi et aDl2008 l). The Log TV-log S re- 



sults from both the OM08 and the current work 
suggest that if the Euclidean nature of the counts 
continue down to the natural confusion as an ex- 
trapolation of our log N - log S, this suggests that 
the cutoff must be fairly sharp below 1 /iJy. 



r.m.s. of ^3.9 ^Jy beam~^ in the best central re- 
gion. It is thus one of the deepest 1.4-GHz radio 
surveys ever undertaken, well matched to the ex- 
cellent Spitzer, Hubble Space Telescope and Chan- 
dra data in GOODS-N, and to the upcoming Her- 
schel, SCUBA-2, and the Large Millimeter Tele- 
scope imaging at 100-850 /Ltm and ~l.l-4mm, re- 
spectively, providing an unbiased probe of star for- 
mation out to z 3 and a means to test the cali- 
bration of other star-formation rate indicators. 

Our analysis of the source size distribution 
and the counts provide further evidence that the 
faint 1.4- GHz emitters have a median angular size 
of ~1.2", that log A^-log S remains flat below 
100 /iJy, and that the natural confusion limit at 
1.4 GHz must be near 1 /^Jy - a prediction that 
can be tested in the coming years using the Ex- 
panded VLA (EVLA) and e-MERLIN. 
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5. Conclusions 

We have presented a deep, new, high-resolution 
radio image of the GOODS-N field and a cata- 
log containing 1,230 secure radio sources in the 
central 40' x 40' region. The image reaches an 
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Fig. 2. — Log iV-log S for co rrected different i al sou rce counts below 1 mJy at 1.4 GHz. HDF B iggs, Lockman 
Hole and ELA IS N2 are fromlBi ggs fc Ivisod (|2006f ) and SSA 13 is from lFomalont et all (|2006l ). The 1046+59 
data are from lOwen &: Morrison (,2008i) . 
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Fig. 3.-- Extended sources in the GOODS-N field. 
Contour levels are listed on the figure and the syn- 
thesized beam (FWHM) is illustrated at the bot- 
tom left. Several of the sources displayed here 
might be extended due to band-width smearing 
issues. The quasi-periodic artifacts are due to un- 
correctable correlator errors. In practice, this pat- 
tern is not the same for every channel and thus 
cannot properly be remove from the images by 
CLEAN deconvolution. 
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Fig. 10. — Brightest source in the HDF is in the 
upper right. Undeconvolvable sidelobes around 
this source are clearly visible and point to the 
phase center. 
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